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Abstract 

Analytical  derivations  are  presented  for  inductor- 
converter  bridge  (ICB)  circuits  in  which  energy  is 
transferred  from  a  storage  inductor  to  a  load 
inductor  with  solid  state  bridges.  These 
derivations  provide  complete  analytical  circuit 
solution  in  contrast  to  previously  available 
numerical  (non-analytical)  procedures.  The 
analysis  is  based  on  two  parallel  methods:  (1) 
Fourier  expansion  of  the  inverter  waveforms  and 
(2)  a  novel  method  based  on  the  inherent  waveforms 
of  the  ICB,  labeled  square  functions.  Our 
analytical  values  of  power  flow,  inductor 
currents,  and  voltages  compare  favorably  with  the 
results  of  a  three-phase  ICB  experiment  at  Argonne 
National  Laboratory. 

Introduction 

The  inductor-converter  bridge  (ICB)  is  a  solid 
state  dc-ac-dc  converter  system  for  reversible 
energy  transfer  between  two  inductors.  This 
system  is  especially  suitable  for  pulsed  power 
supply  applications  greater  than  several  hundred 
megawatts  and  durations  from  a  fraction  of  a 
second  to  many  seconds.  Two  such  applications  are 
the  superconductive  equilibrium  field  coils  of  the 
projected  tokamak  fusion  power  reactors  and 
superconductive  magnets  to  be  used  in  future 
particle  accelerators. 

The  ICB  system  is  inherently  efficient,  control¬ 
lable  in  real  time  and  allows  isolation  of  large 
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pulsed  reactive  loads  from  the  power  grid.  Thus, 
only  the  average  system  losses  are  drawn  from  the 
grid. 

Operation  of  the  ICB  Circuit 
Detailed  operation  of  the  ICB  may  be  found  in 
references  1  and  2.  Figure  1  shows  a  three-phase 
ICB  where  the  storage  and  the  load  coil  is 
represented  by  Lg  and  L^,  respectively.  At  a 
typical  instant  during  the  energy  transfer,  dc 
currents  ig  and  will  be  flowing  in  the  storage 
and  load  coils,  respectively.  The  SCR's  of  the 
left  hand  side  (storage  side)  are  fired  in  the 
normal  Graetz  bridge  sequence:  SL1  SL5 ,  SL1  SL6, 

SL2  SL6»  SL2  SL4>  SL3  SL4>  SL3  SL5>  SL1  SL5’ 

The  SCR's  of  the  right-hand  side  (load  side) 

follow  the  same,  sequence  but  may  be  out  of  step 
with  respect  to  the  storage  side.  The  direction 
and  the  level  of  power  flow  is  determined  by  the 
relative  timing  between  the  source  and  the  load 
bridge  switching  sequences  such  that  a  load  bridge 
lead  will  cause  power  flow  into  the  load  and  vice 
versa.  The  Y-connected  capacitors  in  the  middle 
serve  as  the  intermediate  energy  store  between  the 
storage  and  load  coils  and  they  provide  the 
reverse  voltages  to  commutate  the  inductor 
circuits  from  one  SCR  to  the  next.  Thus,  no 
external  counterpulse  circuit  is  needed. 

Only  a  very  small  fraction  of  each  coil  energy  is 
extracted  in  each  bridge  cycle.  Therefore,  by 
varying  the  relative  timing  between  the  source  and 
load  switching  sequences  and/or  the  frequency  of 
operation,  very  fine  control  over  the  rate  of 
energy  transfer  can  be  achieved.  The  functional 
dependence  of  the  power  on  the  relative  timing 
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(relative  phase  difference)  and  the  frequency  will  The  total  instantaneous  average  power  delivered 
be  illustrated  in  the  following  sections.  from  the  storage  coil  is  m  times  the  above,  or 


Circuit  Analysis  Based  on  Fourier  Components 

Figure  2  shows  an  m-phase  ICB.  For  this  analysis, 

the  SCR's  have  been  replaced  by  ideal  switches. 

This  idealization  implies  a  lossless  system  and 

operation  within  the  successful  commutation 

bounds.  Since  each  bridge  cycle  changes  the  coil 

energies  by  a  very  small  amount,  the  inductor 

currents  are  nearly  constant  in  one  cycle. 

Viewing  such  inductor-converter  half  from  the  m- 

phase  ac  lines,  one  will  see  an  m-phase  square 

wave  current  source  system,  each  phase  being  — 

m 

radians  displaced  from  the  next  and  the  wave 
amplitude  being  equal  to  the  instantaneous  coil 
current,  Fig.  3.  Since  the  average  power  is 
divided  equally  between  the  phases,  it  may  be 
calculated  from  a  one-phase  diagram.  Fig.  4. 

The  calculated  instantaneous  average  power  from 
the  left-hand  source  (storage)  to  the  right-hand 
source  (load)  is 


<PS  (t)>  =Yj 

n 


2m  ISIL 
3  2 

n  it  me 


(-1) 


.  .  nr. 
(sin  — ) 
m 


sin  na 


The  effect  of  the  relative  phase,  a ,  and  the 

bridge  switching  frequency,  <u  ,  on  the  power  flow 

is  evident  in  the  above  relationship.  Note  also 

that  the  contribution  of  higher  harmonics  is 

attenuated  by  the -m-  term.  Figure  6  represents 
nJ 

the  plot  of  the  instantaneous  average  power  as  a 
function  of  ct  with  the  number  of  phases,  m,  as 
the  parameter.  For  m  »  3,  over  99.5%  of  the  net 
power  is  delivered  by  the  fundamental  frequency, 
making  the  m  »  3  curve  in  Fig.  6  almost  purely 
sinusoidal. 

The  time  functions  of  average  coil  currents, 
power,  and  voltages  may  now  be  calculated.  For 
simplicity,  let 


<?Sa  ^■t^>  2nuic  ^anSa  bnLa  bnSa  anLa^ 

n 

where  the  a's  and  b's  are  the  Fourier  coefficients 
of  the  storage  and  load  source  waveforms  as 
indicated  by  the  subscripts.  For  illustration, 
the  Fourier  coefficients  of  the  symmetrical 
waveform  of  Fig.  5  which  is  for  when  m  is  odd  will 
be  used. 
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Then  from 


<p„  (t)>  = 


d  ES(t)  _  d  EL(t) 
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where  a  is  the  angle  by  which  the  load  bridge 
leads  che  storage  bridge.  After  substitution,  the 
average  power  per  phase  will  simplify  to 


Solving  these  equations  with  the  initial 
conditions 


(?ig>  | t=0  =  IQ  =  initial  storage  current 
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we  obtain 
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'<is>  (t)  =  IQ  Cos 


S  L 


-<iL>  M  =  I0a/lJ  S±n  '=  C- 


yjLs\ 


The  average  power  will  become 


<PS  >  (t) 


1/2  K  IQ2 


Figure  8  shows  how  a  phase  current  of  the 
equivalent  three-phase  circuit  may  be  decomposed 
into  two  Sq  functions.  Thus,  the  currents  iga  and 
iLa  of  Fig.  4  may  be  written  as 


o  <  t  <  T,  o  <  T  <  T,  o  <  T  <  1/2 


where  t  is  the  relative  lag  time  of  the  storage 
bridge. 


and  the  average  coil  voltages  from 


<v  > 


<vs>(t) 


<vL>(t)  =  K  I0 


Cos  — ^ 
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Circuit  Analysis  Based  on  Square  Function 
Calculations 

The  preceding  Fourier  method  shows  the  influence 
of  harmonics  in  the  system.  However,  for  real 
time  control,  solving  the  equations  of  the  form 


<P>  =  S  fn  (c0 
n 

for  a  is  time  consuming  and  possibly  inaccurate. 
The  following  technique  will  provide  closed  form 
equations  that  are  efficiently  solved  by  a 
microcomputer  in  real  time  control. 


The  calculations  are  based  on  specially  tailored 
functions  symbolized  by  Sq  (X)  and  Tr  (X),  Fig.  7. 
These  waveforms  are  inherent  in  the  operation  of 
the  ICB  circuits.  The  Sq  function  is  a  good 
mathematical  representation  of  the  ac  phase 
currents  of  the  system.  The  Tr  function  is  the 
integral  form  of  Sq  and  is  a  good  representation 
of  the  capacitor  voltages  due  to  the  phase 
currents.  A  brief  mathematical  development  of  the 
Sq  and  Tr  functions  is  presented  in  the  Appendix. 


These  current  representations  for  one  period  are 
adequate  for  instantaneous  average  power 
calculations.  The  net  power  out  of  storage  per 
phase,  results  from  the  interaction  of  iga  and 
vLa,  the  capacitor  voltage  due  to  the  load  phase 
current: 


if! 
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This  power  is 
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The  instantaneous  average  power  per  phase  is 
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[Sq  <7  +  ?  -  Tq) 


Tr  (t  +  !•)  +  Sq  (x  +  ~  -  t  )  •  Tr  (t)  + 
0  0  0 


Sq  (t  -  T  )  •  Tr  (t  +  t)  +  Sq  (x  -  x  ) 
o  o  o 


•  Tr  (t)]  dT. 

The  four  terms  in  the  integrand  are  first 
transformed  to  the  normalized  variable  functions 
shown  in  the  Appendix,  then  each  term  is  evaluated 
by  using  the  proper  integral  identity  in  the 
Appendix.  The  result  is  then  multiplied  by  three 
for  the  total  three-phase  power 
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where  yq  =  ~  “  jso3  is  *-^e  normalized  lag  time  of 
the  storage  bridge  relative  to  the  load  bridge  and 
the  frequency  is  represented  by  the  period  T. 

This  is  the  closed  form  of  the  <p>  vs  a  curve  for 
m  =  3  in  Fig.  6.  The  time  functions  of  the 
average  coil  currents,  power  and  voltages  may  be 
calculated  as  before.  The  only  difference  being 
that  K  (yq)  is  in  closed  form. 


Expressing  the  actual  circuit  waveforms 
analytically,  allows  other  useful  calculations 
such  as  the  actual  capacitor  voltages  throughout 
the  transfer  cycle,  the  study  of  commutation 
throughout  the  cycle,  and  the  actual  coil  voltages 
and  instantaneous  currents,  without  resorting  to 
numerical  procedures. 


Comparison  with  Test  Results 

A  model  three-phase  XCB  has  been  built  and  tested 
at  Argonne  National  Laboratory.  This  system  uses 
two  identical  superconducting  coils  capable  of 
storing  125  kj  at  250  A,  as  the  storage  and  load 
inductors.  Other  system  parameters  are: 


Ls  =  Ll  =  4  H,  IQ  =  iOO  A 
C  =  10-4  F,  a  =>  90° 


<v^>  (t)  =  224.4  cos  0.5609  t 
<ps>  (t)  =  11218  sin  1.1218  t 

A  plot  of  these  equations  appears  in  Fig.  9. 

Figure  10  shows  the  average  coil  voltages  and 
currents  obtained  experimentally.  Good  agreement 
exists  between  the  analytical  and  the  experimental 
results.  Note,  however,  that  the  experimental 
initial  storage  current  is  somewhat  higher  than 
the  final  load  current.  This  is  due  to  the  losses 
in  the  system  which  is  neglected  in  this  analysis, 
but  may  be  incorporated  in  the  differential 
equations  leading  to  the  average  time  functions. 

Conclusions 

The  behavior  of  the  multiphase  inductor-converter 
bridges  have  been  studied  by  two  analytical 
techniques.  The  conventional  Fourier  technique 
produces  the  average  circuit  power,  currents,  and 
voltages  as  a  function  of  time.  It  also  shows  the 
effect  of  the  existing  harmonics  in  the  circuit 
behavior.  The  square  function  technique  is 
particularly  devised  for  the  ICB  and  other  SCR 
circuits  in  which  rectangular  waveforms  appear. 

The  identities  defined  on  the  special  functions 
Sq  (X)  and  Tr  (X)  operate  directly  on  the  circuit 
waveforms.  Thus,  much  more  information  about 
instantaneous  behavior  of  the  circuit  is  available 
for  analysis,  system  design,  and  the  development 
of  real  time  control  algorithms.  Preliminary 
tests  with  open  loop  microcomputer  control  have 
been  conducted  with  satisfactory  results.  The 
development  of  an  optimal  closed  loop  control 
algorithm  is  currently  in  progress. 


u  =  4084  rad/s 

The  system  equations  are  derived  from  substitution 
of  these  parameters  and  m  =  3  into  the  time 
functions  shown  in  the  Fourier  analysis  section: 

<ig>  (t)  =  100  cos  0.5609  t  A 

<i^>  (t)  =  100  sin  0.5609  t  A 

<vg>  (t)  =  224.4  sin  0.5609  t  V 


Appendix 

The  Sq  function  is  defined  as  the  sum  of  unit  step 
functions  as  follows, 

Sq  (Y+Y0)  £  u  (y)  -  2  u  (y  -  y  +  Yq)  +  2 
u (Y  “  1  +  YQ)  -  u  (yq  -  1)  , 

0  If  1  I  ,  0  <  Y  i  1, 

O 
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The  Tr  function  is  defined  as  the  integral  value 
of  the  Sq  functions, 

Tr  (y+Tq)  &  f  Sq  (c  +  y  )  d? 

•J  o 


Tr  (t  +  t  )  -  f  ^ 

0  Jo 

where  Y  =  £  and  Y„ 
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_  t£ 

T  ° 


0  <  t  <  T 


0  t  <  i 

o  2 


“  Y  u  (y)  -  2  (y  -  -r  +  v  ) 

Z  O 

U  (Y  “  W  +  Y  )  +  2  (y  —  1  +  Y  ) 

C*  0 

u  (y  -  1  +  Yq)  -  (y  -  1)  u  (y  -  1) 

The  following  integral  is  extensively  used  in 
calculating  the  average  power  flow  in  the  circuits 
of  interest.  Therefore,  it  will  be  stated  as  an 
identity  which  may  be  directly  verified. 
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J  n 
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[Sq  (y  +  Yx) 


Tr  (y  +  Y,)]  dY 
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2y2  ~  Y1  “  4y1y2  +  y2 


v  +  2y2,  yx  >  y2  > 
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The  following  identities  will  also  be  of 
considerable  value 


Scl  (Y  "  Y0)  “  ~  Sq  [y  +(l/2-y0)], 

Tr  ( y  —  y0)  -  -  Tr  [y  +  (l/2-y0)  I ,  o  <  yo  <  i 

which  may  be  verified  by  direct  substitution. 

Note  that  by  using  these  identities,  we  can  easily 
evaluate  the  above  integral  for  any  combination  of 
leading  and  lagging  functions. 

The  above  Sq  and  Tr  functions  have  been  defined 
for  a  period  equal  to  1.  For  periods  other  than 
1 ,  the  argument  should  be  multiplied  by  the 
appropriate  constant: 


Sq  (t  +  tQ) 


Sq  (y  +  v  ) 
■  o 


{0  <  t  <  T,  0  <  t  <  £ 
°  2 

0<y<l,  0<yo<i 


Substitution  of  y  and  yQ  in  the  above  integral 
will  give 


Tr  (t+tQ)  =  tJY  Sq  (5  +  YQ) 


dC 


=  T  [Tr  (y  +  yQ)]  =  t  u  (t)  - 

(t  -h  v  v  (t  -f  +  to) 


+  2  (t  -  T  +  t  )  u  (t  -  T  +  t  ) 
o  o 


(t  -  T)  u  (t  -  T). 


For  time  base  calculations,  these  functions  may  be 
used. 
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Fig.  1.  Circuit  Diagram  for  the  3- 
Capacitor  Model  IC  Bridge. 
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Fig.  2.  Circuit  Diagram  for  an 
a-phase  XCB. 
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Fig.  3.  Equivalent  Diagram  of 
an  m-phase  ICB. 


Fig.  4.  One  Phase  Equivalent 

rircuit  of  an  m-phase  ICB. 


Fig.  5.  Single  Phase  Current 
Source  Waveform  for 
an  m-phase  ICB  Circuit. 
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(a)  Phase  Current  in  a 
Three-Phase  ICB 

i  (b)  Decomposition  of 
the  Phase  Current 
into  Two  Sq 
Functions . 
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Fig.  9. 

.  Average  Coil  Currents 
I-  Voltages,  and  Power 
vs  Time. 
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Fig.  10. 

Experimental  Average  Coil 
Voltage  and  Current  Traces. 
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